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Abstract 



We consider the effect of a non-zero lattice spacing on the low-energy effective theory of Wilson fermions with Nf = 1. Analyt- 
ical results are given for both the chiral condensate and the microscopic spectral density of the Wilson Dirac operator. A subtle 
' r _ m ^ "mechanism ensures that a constant chiral condensate is recovered, once the sum over sectors of fixed index v is performed. 
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1. Introduction 

The low energy effective theory of QCD is based on Gold- 
stone bosons of the broken chiral symmetry S U(Nf)L ® 
SU(Nf)R. As is well known, the t/(l) axial symmetry is ex- 
— i plicitly broken at the quantum level, because of the anomaly. 
3 This means that the Nf — 1 theory might appear quite trivial at 
' low energy because there are no Goldstone bosons. Leutwyler 
^"'and Smilga showed in [Q]] that this is not the case: an expansion 
of the partition function in terms of me' 6 for Nf = 1, shows that 
1 interesting behavior arises in sectors of fixed topology v. 

For simulations of QCD on the lattice, one faces the problem 
of doubling. For one flavor, the 15 extra states can be removed, 
by adding a Wilson term to the Dirac operator. We denote Dirac 
operator with the Wilson term included by Dw- The addition of 
the Wilson term, moves the eigenvalues of the extra states away 
from the origin in bunches of 4, 6, 4, and 1. In continuum 
language the Wilson term is equivalent to a double derivative 
proportional to the lattice spacing a. 

The Wilson term destroys the anti-hermiticity of the Dirac oper- 
ator, only retaining y 5 -helicity, D f w - ~y 5 Dw~y 5 - The hermitian 
operator D$ - ~y 5 (Dw + m) can still be defined. The eigenval- 
ues for Dw spread in the complex plane, removing the apparent 
definition of a sector v in terms of zero modes. 

We explore the low energy behavior of the Wilson Dirac 
operator by using the low-energy effective field theory, for 
which we confine ourself to the e regime. For the Nf = 1 
case, there is no spontaneous breaking of symmetry, but it is 
still possible to find a behavior very closely reminiscent of the 
e-regime in chiral perturbation theory, if we focus on sectors 
of fixed index v. A choice of definition of such sectors as the 
fourier modes of the partition function, ensures a simple defi- 
nition of v, as the average sign of chirality of the real modes J^] . 

We will explore the special case of Nf = 1. We start in 
section 2 by remembering the effective theory for low-energy 
Dirac operator for Nf = 1, and in section 3 we include the 
terms emerging from the Wilson term. In section 4 we show 
how to obtain the chiral condensate, which will be used to find 



X 



both the full condensate and the condensate for each sector 
of fixed index. We use this to demonstrate a quite intriguing 
behavior when the lattice spacing is non-zero. We will in 
particular focus on how the full condensate is built up from 
sectors of fixed v. In section 5, we look at the spectral density 
of the hermitian Wilson Dirac operator D$, as well as the 
corresponding p x , where we include two terms not previously 
considered in the Nf — 1 case. 

2. Nf in the continuum 

Here we recall the arguments of Leutwyler and Smilga 1 1 ] 
for Nf - I. It was shown, that though there are no Gold- 
stone bosons for Nf - 1, the effective partition function in the 
e regime, where one confines to a box of size L - - and 
l^, could still be described by the exponential to minus the en- 
ergy 



Z(m, 9) 



YVmRe(e M ) 



(1) 



By use of the U(l) axial anomaly, it was shown that the only 
parameter the energy would depend on, was me' 6 , for which the 
given result is the lowest order in m. The decomposition of a 
sector of fixed index in QCD 



Z v (m) = (m M U k (A 2 k + mm*)), 



(2) 



was used to find the transformation properties. At are the eigen- 
values of the Dirac operator. These transformation properties 
could be mimicked in the effective partition function if the 
fourier modes of Z, was chosen as the effective partition func- 
tion for fixed index v, ie. Z v such that 

Z(m,8) = ^V flv Z>i). (3) 



3. Inclusion of Lattice Spacing 

We now extend to non-zero a. The e counting scheme for 
a + 0, is taken to be 



v- 



(4) 
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such that first order terms in m are compatible with second or- 
der of a. The terms proportional to the lattice constant a are 
included in the effective theory, by the same principle as how 
the mass is included, as they both comes from terms depending 
on <A*A 0]- 

VL a = Vc Q Tr(a(U + C/ t )) + VW 6 Tr(a(U + t/ 1 ")) 2 (5) 
+ VW 7 Tr(a(U - t/ 1 ")) 2 + VW s Tr(aUaU + af/W + ). 

The first term is of the same form as the first order term in m 
and thereby redefine both terms to ^Tr(M(U + t/' )), where we 
have defined VLm + VW&a = M. The remaining terms are set 
to a 2 W 6 V = a 2 , a 2 W 7 V = a 2 , a 2 W & V = a 2 , which are all ~ 1 
in the e regime. For a general Nf with Goldstone bosons, this 
gives |H] 

dUei Trmu+ui)) 

iU{N f ) 

Xe -alTr(u+u i ) 2 -a?jTr(u-ir) 2 -alTr(u 2 +(u f ) 2 ) 



(6) 



For Nf = 1 with no Goldstone bosons, we obtain the same form 
by adding terms with ae' 6 , to the partition function, such that we 
only obtain real terms. We have chosen to add the phase to the 
lattice spacing a. In the e regime this gives 

Z(8,m,a) = exp(racos(fJ) - 2a 2 ; cos(2f7)). (7) 

The parts proportional to a 2 6 and a 2 have been omitted, because 
for Nf = 1 for the partition function, a?, a 2 and cp- % can be re- 
defined into a a 2 , term. If the fourier modes of Z are chosen as 
a definition of the index on the lattice, then according to a 
sector v will be given by 



(8) 



which is only effected by the real eigenvalues A, since A with 
Im[A] + then (k\y 5 \k) = 0. With the inclusion of the 1/(1) 
integral, Z v for Nf = 1 becomes 



f' 



Z v (m, a 8 ) = I d8exp{iv0 + m cos(6») - 2a\ cos(2#)). (9) 

r 

This defines the 9 and v for a + which we will use. 

4. The Chiral condensate for Nf = 1 

We are interested in the analytic behavior for one flavor Nf = 
1 with a + 0. The Chiral condensate is found from the partition 
function 



2 8 2 

£(m, a%,0) — — — Z(6, m, a 8 ). 
For one flavor using (7) we obtain 



a s , 6) 



cos(0), 



(10) 



(ID 



which is exactly as found in [6] in the continuum. This result is 
in contrast to the much more complicated condensate for each 
sector v, 57(m, a%), which is given by 



2 v (m,«8) = 



j"^ ddcos(ff) exp(iv0 + m cos(0) - 2a^ cos(26')) 
j_ d6 exp(iv8 + m cos(ff) - 2a^ cos(20)) 




Figure 1 : The mass dependent Chiral Condensates Z„ for Nf = 1. The graphs 
shows the results for v = 1 and v = 2 for «g = 1.5. I is observed how the 
additional effective term from the non-zero lattice spacing, pushes the v poles 
away from the origin and the v — 1 zeros for |v| > 2 and out onto the real axis. 



Our first observation concerns the divergences of E v . Making 
the change — » -9 in both integrals, it is seen that Z v = £~ v , 
and we are therefore only interested in the behavior for positive 
v. When a — 0, there is a divergence at m = for v + 
coming from the zero-modes. When a non-zero lattice spacing 
is included, this gets smeared out as in figureQ] The divergence 
at the origin is now only found for odd v. This is dictated by the 
symmetries of Z v since 



Z v (m, a%) 



- f 

Jo 

- f 

J -71 



givG gtncos(6)-2a% cos(2f9) 



Q inv giv6 ^-m cos(f?)-2flg cost 



(13) 



= Z"(-m, fl8 K-l) v . 

The remaining divergences are pushed further out on the real 
axis, where they appear in pairs symmetric around m — 0, 
which is seen from the symmetry used to find the divergences 
for the odd term. Expanding in «g << 1 and m < 1, one finds 
the behavior to be 



? v=2 



r v=3 



(m/2 - a\ + w 3 /12) 
m 2 /4 - 2a 2 - a 2 m 2 /2 + m 4 /48 
1 -a 2 + m 2 /8 -a 2 m 2 3/8 
m -a\ + m 3 /24 - ajm 3 / 16' 



(14) 



(12) 



where the pole for the odd term in m = can be seen. The poles 
position for small m are apparent from this approximation, for 
which the first pole are at m oc a 8 . The behavior of the zeros of 
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Z v are plotted in figure[2] We see that the real poles have a linear 
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Figure 2: Divergence plot for 2 V , The plot shows a 3d plot for Re(Log(Z v )) 
seen from above for v = 3 and 7 respectively, where the jc-axis is ajj and the 
y-axis is m. The divergences are controlled by the zeros of Z v and the real part 
of the logarithm blows up at Re(Z y ) = Im(Z y ) = which is seen as the lines in 
the plot. 

dependence for small a?, while for bigger a% the dependence 
becomes quadratic. 

We are also interested in finding how each sector add to the 
full chiral condensate 



ln( j>] e ive Z v ) 



(15) 



V— — CO 

„iv8' 



d_ 

dm 

Lav=-co c u flI A 

yv=co i v e 7 v 

ZjV=-O0 e ^ 

CO 

V e iv6 I. v Z v /Z. 



This tells us that we should weight each sector v with the factor 
e' vfl Z v /Z, such that Z v cancels the poles. In figure[3]we see how 
for 6 — this sums up from -v to v. The higher terms cancel 
the too high values around m = 0, and this cancellation is too 
big, such that X at m = oscillates, though converging towards 




Figure 3: jj y= Z V V /Z for a\ = 1.5 and ; = 0, 1,3,9. blue(J = 0), Purple(j = 
1), Yellow(j = 3) and Green(j = 9). The oscillating behavior for «g # is 
observed. 



w 0.6 




Figure 4: £:£ = _ . Z V Z"/Z for = and = 0, 1,3, 9. blue(j = 0), Purple(j = 
1), YellowQ = 3) and Green(j = 9). We compare how the sum over fixed 
sectors built the full condensate. 



X = 1. We compare this to the a — case in figure |4] We 
find that the sector v = 0, is not too significant, as there is no 
contribution to X at m = 0, and X° in no way fills up most of X. 
That X° is almost 1 at m — 1 in figure[3]is a coincidence and is 
not true for other a?. 

While we see how different the picture is, it is important to note 
that this is mainly if one focus on the pole behavior around 
m — 0. As one passes the poles at around m = 8flg the behavior 
quickly goes to that of ag = 0. We show this in figure|5] 

5. The density with a* and a* 



1, we redefined a 2 6 and a 2 into a?, such that Z v was 

2 



For N/ 

only dependent on at. For the spectral densities of the Wilson 
Dirac operator, this is not the case. As in |2|], we expand the 
partition function to 



Z m (M,Z,ad 



det(D w + nif + z/y 5 )det(Diy + m + zys)\ 
det(LV + m' + z'ys) j 



where we have included the ys -mass term, which is propor- 
tional to Z. We have added one extra bosonic flavor and one 



3 



which was done without and a 2 in (§]. The terms with and 
a 2 are a bit more complicated, but can still be reproduced here 

(19) 



S f — cos(6)(-nif sin(f + u) - m sin(f - u) 
+ iz/cos(t + u) + 4flg cos(2f) cos(2m) cos(#) 

2 
8 

6 



+ fe cos(f - u)) - Aa\ cos(2f) sin 2 (f/) 
S b — -im' sinh(s) - iz' cosh(s) - 2a\ cosh(2s) 
$67 = 4a?(2/ cos(rJ) cos(m) sin(f) - sinh(i)) 2 



10 



15 



2(1 



Figure 5: £ (as = 0)/£-(cig = 1). The ratio for the Chiral condensate between 
the continuum and a lattice effective theory for a fixed sector v = 2, is seen to 
converge to 1 , when m is away from area of the poles around the origin. 



more fermionic flavor. The partially quenched condensate is 
then defined as 

S y (m/,m) = 3,„(ln(Z ( v 2|1) (m/,m,m'))| m-m ■ (16) 

It was shown in ^\ that spectral densities could be obtained by 
taking the discontinuity, which for zZ v {m /, m) obey 



2>, 



" 1 A w )sign[(k\y 5 \k)] 



(17) 



+ 4a 2 (2 cos(fj) cos(m) cos(f)) - cosh(i)) 2 

Pn = i— (cos(0)(z'm/sin(£ + u) + z/cos(t + u)) + m sinh(s) 

+ z' cosh(s)) + a\(2 cos(2f + 2k) cos 2 (6>) - 2 cos(2f) sin 2 (6>) 
+ 4 cosh(;f + iu + s) cos(fT) + 2 cosh(2s)) 

+ 4«g(2/ cos(fT) cos(m) sin(f) - sinh(s))(/ sin(f + u) cos(ff) + sinh(s)) 
+ 4a 2 (2 cos(ff) cos(m) cos(f) - cosh(s))(cos(? + u) cos(fT) + cosh(s)) 

P22 - i-(cos(9)(im sin(f - u) + zcos(f - u)) + m sinh(s) 

+ z' cosh(s)) + «g(2 cos(2f - 2u) cos 2 (0) - 2 cos(2f) sin 2 (6>) 
+ 4 cosh(;f - iu + s) cos(fT) + 2 cosh(2s)) 

+ 4«g(2/ cos(6) cos(m) sin(f) - sinh(i))(/ sin(f - u) cos(ff) + sinh(s)) 
+ 4a 2 (2 cos(f?) cos(m) cos(f) - cosh(s))(cos(f - u) cos(fT) + cosh(s)) 



-- sin(6*)(m/e 
4 



■ z f e-" + ze") 



= -Im[2Z v (m f ,m)]\ m=m , =A w, 
n 

where sign[(k\ys\k)] were included for regularization. /lF are 
the eigenvalues of D\y- The extra quarks result in an en- 
larged symmetry group, and we therefore need to integrate 
over the largest convergent sub-group of G/(2|l) as described 
in Q. While the Grassmann integration trivially converge, the 
bosonic integrals need careful attention. Following [2] we ro- 
tate U — > ill to make Z v convergent for positive a?, such that 



- 4«g sin(f7)(sin(? - is) + sin(2f) cos(m) cos(fT)) 

+ 4flg(2/ cos(fT) cos(m) sin(f) - sinh(s))/cos(f) sin(fJ) 

- 4a 2 (2 cos(6>) cos(m) cos(f) - cosh(s)) sin(f) sin(fT) 
1 



• sin(0)(m/e B + me " + Zf e " ~ ze ") 



"(21 1) 



-K2 



Z£ m = I dUSdet(U) v exp(-Str(M(U - CT 1 )) + a\Str(U - IT 1 ) 



- 4ag sin(f7)(sin(f - is) + sin(2f) cos(m) cos(fT)) 
+ 4flg(2/cos(6') cos(m) sin(f) - sinh(i))/cos(f) sin(f?) 
4a 2 (2 cos(6>) cos(m) cos(f) - cosh(i)) sin(f) sin(fJ) 



+ a$Str(U+ U~ l ) 1 + a\Str{U L + U~ z ) + -Str(Z(U + U~ 1 ))). 
We choose the parametrization 





' e i,+i " cos(6») 


ie it+i * sin(0) 





s 




' 







u = 


ie"- ' 1 ^ sin(0) 


e u ~ iu cos(0) 







exp 








C2 




I 





e s 






I p\ 




J 



P4 = — cos(6)(imf sin(f + m) + im sin(f - m) + zy cos(f + m) 

+ z cos(f - m)) h (m' sinh(s) + z' cosh(s)) 

1 2 
+ flg(- cos(2f) cos(2m) cos 2 (0) + cos(2f) cos 2 (0) + - cos(2f) sin 2 (#) 

8 4 
+ - cosh(/f + s) cos(m) cos(0) + — cosh(2i)) 



for which the Berezinian becomes J = 4cos(#) sin(0)(l + 
\(aiPi + arfi)) 0] and we set |7| = 4| cos(6») sin(6»)|. Plugging 
U into the partition function Z v we find that the <p dependence 
disappear, and the <p integral simply gives n. The Grassmann 
integrals can also be carried out explicitly. This gives 



2 1 

+ 2ag(-(2/cos(f?)cos(M)sin(f) - sinh(s)) 



Z 2|l( m /' m > m '>Zf> 2, z, flg) 

Xoo y-vr r*n r*n 

OO xj -71 %J—7T *J —IX 



(18) 



S f +S b +S 6 7+(2it-s)v 



X (P 4 - ^11^22 + ^12^21 - -(Pit + PllWl 



x (/ cos(ff) cos(m) sin(f) + sinn^?)) 

- (/ sin(f + u) cos(ff) + sinh(s))(/ sin(f - u) 
x cos(f/) + sinh(i)) - sin(0) 2 cos(f) 2 ) 

? 1 

+ 2a^j(—(2i cos(fT) cos(w) cos(f) - cosh(i)) 

x (cos(fT) cos(m) cos(f) + cosh(s)) 

- (cos(f + u) cos(fT) + cosh(i)) 

x (cos(f - u) cos(0) + cosh(s)) + sin(fj) 2 sin(f) 2 ), 



for which convergence requires > a b + - a*, 
representation can be found in We numerically solve the 
rest of the integrals to obtain the densities. In figure [6] we show 
a couple of densities for Nf = 1. It is seen how the addition of 




Figure 6: p v x (A 
Purple(a 2 6 = 0.01, a* 
and Green(a^ = 0, 



,nif) for v = 2 and irif = -1. Yellow(a^ = aij = 0, «g = 0.04), 
a\ = 0.04), B/He(«2 = a 2 = 01fl 2 = Q4) 

flg = 0, alj = 0, c/g = 0.05). The density p" x for specific sectors v 
compared for different terms in the low-energy effective theory. 



the and a?j term changes the positions of the tops and changes 
the values. It should be noticed how the density vanishes at 
m = — 1, as is required from the QCD partition function since 
p v x is an average of J^ k 6(A W + A^)Ylj(Aj + nif), which is always 
zero for nif = m - A w . 

When we defined the sectors v for the partition function, we 
said that a sector v is equal to the average of chirality for a 
configuration. For the density p v x this is equal to 



r Px (A w )dA w 

U -co 



6(Af + A w )sign[(k\y 5 \k)}) dA 



J-oo \ k 

= {^ J ^8n[(kMk)Vj =v. 



(20) 



As a central self-consistency check, it has been verified by nu- 
merical integration, that this is true for «g, and a\ at sizes like 
in figure [6] 

Finally we also consider the spectral density of the hermitian 
Wilson operator D5 = ys(Dw + m) obtained as in Q by 

p v 5 (A 5 ,m f , a) = 6(A 5 k - A 5 )^ = ^Im[G v (-A 5 ,m f , a)], 



(21) 



G v (z,m f ,a) = ( — \nZ* m (M,Z,a)\\ z =z' ,m f =m=m> ■ (22) 



where G(z) is the resolvent of D5 obtained as 

d_ 

dz 



z and z' are the 75 mass which we included in the graded parti- 
tion function, which comes from the term ztfryst//. Zf is always 
set to 0. The result of including the and term, can be seen 
in figure Q and [8] We see that p v 5 is at the origin as required, 



1.2 



i.n 



0.8 - 



0.6 - 



0.4 



0.2 



0.0 









/ v = : 









,15 
= 1. 



0,aj 



0.04), 



0.01, a\ =0.04) and 



Figure 7: p^(A ) for v = and nif = 1. Blue(a£ = 
Yellow(a\ = 0.01, azj = 0,a\ = 0.04), Green(a\ = 0,a?j 
Purple(c? 6 = Q,a?j = 0, ajj = 0.05). The Density p* for fixed sectors v compared 
for different terms in the low-energy effective theory. 
1.4 




Figure 8: p v s (A ) for v = 1 and my = 1. Blue(a 6 
Yellow(a 2 = 0.02, a] = 0, ajj = 0.05) and Purple(al 



a] 



= 0,4 
0.02, al 



0.05), 

: 0.05). 



The Density pi for fixed sectors v compared for different terms in the low- 
energy effective theory. 



since p v 5 depends on the average of Y*k8(A 5 - A 5 k )Ylj(A 5 j) which 
will always be at the origin. 

Very recently, two lattice QCD studies have demon- 

strated that the quenched predictions for the microscopic spec- 
tra from WCPT can be matched to the lattice data. 



6. Conclusion 

We have considered the effect of a non-zero lattice spacing, 
in the low-energy effective QCD partition function, for one fla- 
vor in the e regime. With this we explored the behavior of the 
chiral condensate for QCD. We saw how this meant that for 
Z v the pole at the origin, coming from the v zero eigenvalues, 
were spread out into v poles. We also saw how the full conden- 
sate was build up. This we compared to the continuous case, 
and saw how the full condensate, for non-zero lattice spacing, 
started to make a damped oscillation, when we summed from 
-V to v. We ended with showing the e regime spectral densi- 
ties of the Wilson Dirac operator, and especially focused on the 



5 



effect of flg and cSj. We saw how a? and had a very simi- 
lar effect of shifting the height of the peaks. We have checked 
that the densities remains zero at the points dictated by QCD. It 
was also checked that f p*(A^)dA v = v, this shows that the 
choice of Z v is consistent with it as the index of Dw- 
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